Regge calculus: a unique tool for numerical relativity 
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The application of Regge calculus, a lattice formulation of general relativity, is reviewed in the 
context of numerical relativity. Particular emphasis is placed on problems of current computational 
interest, and the strengths and weaknesses of the lattice approach are highlighted. Several new and 
illustrative applications are presented, including initial data for the head on collision of two black 
holes, and the time evolution of vacuum axisymmetric Brill waves. 
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I. NUMERICAL RELATIVITY 

The complexity of the Einstein equations, combined 
with the sparsity of relevant analytic solutions, necessi- 
tates a range of other tools with which to explore com- 
plex physical scenarios. These include series expansions 
and perturbation techniques, together with numerical so- 
lutions of the fully non-linear equations. Unfortunately, 
the numerical solution of Einstein's equations has proved 
to be an exceedingly difficult problem. Over the last 
three decades numerous schemes have been developed, 
or adapted, to tackle the vast range of problems which 
fall within the purview of numerical relativity. 

The classic approach, involving a three-plus-one di- 
mensional split of space and time, was first expounded by 
Arnowitt, Deser and Misner (ADM) "T]. This approach is 
natural in the context of our Newtonian intuition, and is 
also directly applicable to computer simulations. Whilst 
the traditional ADM approach has dominated numeri- 
cal relativity, questions regarding its long-term stability 
have lead to the development of many other formulations 
of the Einstein equations in three-plus-one dimensions. 

Techniques of current interest include those developed 
by Sasaki and Nakamura, and later by Baumgarte and 
Shapiro, known generically as Conformal ADM (CADM) 
formulations. Incorporating insights from York's analysis 
of the initial value problem, these algorithms have shown 
superior stability properties compared with ADM in 
some applications. Other recent formulations are based 
on symmetric-hyperbolic forms of the Einstein equations, 
where it is hoped that the mathematical proofs of stabil- 
ity and well-posedness confer numerical advantages com- 
pared with techniques whose mathematical structures are 
more uncertain. The recent review by Lehner Q dis- 
cusses many of these issues. 

Despite the enormous effort invested in these tech- 
niques (and many others, including characteristic for- 
mulations), most modern numerical relativity codes con- 
tinue to suffer from problems with long-term stability 
and lack of accuracy. Lack of resolution is only a partial 
solution; new insights or techniques seem to be required 
to overcome many of these problems. 



II. REGGE'S FORMULATION OF GENERAL 
RELATIVITY ON A LATTICE 

Regge calculus is a formulation of general relativity 
on a piecewise flat simplicial complex, rather than a dif- 
ferentiable manifold. In general this complex is built from 
four-simplices, the four dimensional generalisation of tri- 
angles and tetrahedra. The interior of each three and 
four dimensional lattice element is intrinsically flat, with 
curvature concentrated on the two-dimensional faces (tri- 
angles). In n-dimensions, curvature resides on lattice el- 
ements of co-dimension two. 

Regge calculus would appear to be ideally suited to 
numerical simulations; it is an inherently discrete formu- 
lation of general relativity, and complex topologies are 
easily incorporated. In addition to its facility in classical 
numerical relativity, Regge calculus has also seen appli- 
cation in the sum-over-histories formulation of quantum 
gravity. In this review we concentrate on the application 
of Regge calculus to classical numerical relativity^etails 
on quantum applications can be found elsewhere [j| . The 
remainder of this section is devoted to the description of 
classical gravity on a simplicial lattice. 

Given a lattice spacetime, which consists of a list of 
vertices, the corresponding connectivity matrix and the 
lengths of all edges, the Gaussian curvature k can be 
locally defined by the parallel transportation of a test 
vector about a closed loop. Calculating the angle through 
which the vector rotates, we obtain 
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The first equality is the definition of Gaussian curvature; 
the second provides the equivalent expression on the lat- 
tice. The natural loop to choose, as indicated in fig- 
ure ^^a), is defined by the area dual to the triangle on 
which the curvature resides. This is a portion of the dual- 
lattice, and may in principle be calculated for a given tri- 
angulation. The test vector is rotated through an angle 
e in the plane orthogonal to the triangle (that is, in the 
plane of the dual area A*). We refer to e as the "deficit 
angle" , the lattice representation of the curvature con- 
centrated on the triangle. The deficit angle is computed, 
given the edge lengths of the lattice, using 
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which is the lattice equivalent of the Palatini identity 
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FIG. 1: (a) A triangle 012 embedded in a simplicial complex 
(not shown). To examine the curvature about 012 one trans- 
ports a vector through a closed loop in the plane orthogonal 
to 012. The simplicial complex itself naturally defines this 
loop - the dual triangle, denoted 012*. (b) An illustration of 
the hyper-dihedral angle 0oi2, measured between the tetrahe- 
dra 0123 and 0124 which hinge on 012 within the four-simplex 
01234. 



where the summation is over all four-simplices k which 
contain the triangle, and 9k is the hyper-dihedral angle 
between the two tetrahedral faces which hinge on the 
triangle within k. Figure ^b) shows this structure for 
the simplex 01234 which hinges on the triangle 012. The 
form given above applies for spacetimes with Euclidean 
signature; similar expressions apply for a spacetime with 
signature — h ++. 

Regge derived the simplicial equivalent of the Hilbert 
action [3| for a lattice spacetime. 
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where the summation is over all triangles i, is the area 
and Ci is the deficit angle of the ith triangle. This sim- 
plicial expression can be intuitively understood by not- 
ing that on a simplicial complex the Hilbert integrand 
has compact support on triangles, each triangle has the 
volume element dV = \/~g d'^x w AA* and the scalar 
curvature is i? « k = e/A*. A more thorough proof in 
this spirit has been provided by Miller 0. 

Lattice edges are the discrete representation of the 
metric, and thus any geometric quantity may be calcu- 
lated once all edges are known. In direct analogy to the 
continuum, the equations of motion are obtained by a 
variational principle; the action is varied with respect 
to the independent variables. This yields the vacuum 
"Regge equations" 
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where Ir is the lattice form of the Hilbert action. These 
equations of motion are implied by the Regge identity 
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From the derivation of the lattice equations it is clear 
that there is a single equation for each edge in the lattice. 
In general we expect at least seven edges per vertex in 
a three-dimensional simplicial lattice (triangulating R^; 
other topologies may change this average slightly) com- 
pared with the six independent components of the three- 
metric per point in the continuum. This can be under- 
stood by comparison with a tetrad formulation, where 
one introduces an extraneous variable and a correspond- 
ing consistency condition. We expect that the lattice 
equations themselves provide the additional conditions. 

From this basic structure it is possible, in principle, 
to construct the simplicial counterparts of any geometric 
object of interest: Kfj_i,, Raf3^iv and so forth. However, 
careful averaging is required if these simplicial definitions 
are to converge smoothly (and pointwise) to their contin- 
uum counterparts. 

Finally, we note that the theory of Regge outlined 
above is one among many conceivable approaches to "lat- 
tice gravity". A related method has recently been de- 
veloped by Brewin 7], where the metric is locally con- 
structed in Riemann normal coordinates from an under- 
lying lattice structure. Brewin is then able to use any of 
the standard formulations (ADM, CADM, etc) to evolve 
the lattice. These approaches are related in general struc- 
ture to the finite element and finite volume methods in 
wide use in computational fluid dynamics and engineer- 
ing. It is not yet entirely clear how lattice approaches to 
general relativity relate to these well-developed numeri- 
cal techniques, but such an understanding would provide 
a useful bridge between the lattice theory and continuum 
methods. 



III. (3+l)-DIMENSIONAL LATTICE GRAVITY 

The (3 -f l)-dimensional formulation of Regge calcu- 
lus is reasonably well understood, although to date there 
have been few applications. Significant new studies will 
be required to obtain a full understanding of such issues 
as convergence and long-term stability. 

The initial value problem is a vital precursor to the 
generic evolution problem. A general technique for solv- 
ing the initial value problem in the thin-sandwich, (3+1)- 
dimensional formulation of Regge calculus has been de- 
scribed by Gentle and Miller 8] . Their approach is based 
on an identification of the geometric degrees of freedom 
and conformal structures employed in the York initial 
value formalism 9] . By associating lattice elements with 
the continuum conformal structure. Gentle and Miller de- 
scribe one possible (although far from unique) approach 
to the construction of two-slice simplicial initial data. 
They successfully benchmark their approach on the Kas- 
ner cosmology Q. 
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Once initial data has been constructed a consistent 
evolution scheme is required. The Regge equations, a set 
of coupled non-linear algebraic equations, can in princi- 
ple be evolved by iteratively solving the entire coupled 
system at each timestep. However, it was recently re- 
alised that the lattice can be structured in such a way as 
to allow the implicit, parallel, decoupled evolution of sets 
of vertices JJlJ. This algorithm is known as the "Sorkin 
evolution scheme" . 

A general (3 -I- l)-dimensional lattice may be con- 
structed from a given three-dimensional lattice by the 
"evolution" (not true geometrodynamic evolution) of in- 
dividual vertices. Each vertex is "dragged forward", off 
the initial hypersurface, and in the process addition edges 
are created to join the "evolved" vertex to its counterpart 
on the initial hypersurface. This process is repeated until 
all vertices have been "evolved" . In this way the connec- 
tivity of the initial hypersurface is replicated, whilst the 
lattice structure between the two surfaces is constructed 
naturally as part of the algorithm. A lattice of this type, 
referred to as a "Sorkin triangulation" allows the 

parallel evolution of individual vertices described above. 

General relativity is fundamentally coordinate invari- 
ant, although one must choose a convenient frame in 
which to perform numerical calculations. In the stan- 
dard (3 -I- l)-dimensional ADM formulation, this is en- 
coded in the freedom to lay down coordinates on the ini- 
tial Cauchy surface, together with the freedom to choose 
how these coordinates are propagated to each future time 
slice (lapse and shift freedom). In a lattice formulation 
these correspond to the freedom to distribute vertices on 
the initial hypersurface and the freedom to propagate 
those vertices, respectively. The latter freedom requires 
the provision of four conditions per vertex during each 
evolution step to determine the propagation of the ver- 
tices; the simplicial lapse and shift freedom. 

Regge calculus has been successfully applied in (3-1-1)- 
dimensions in the context of the anisotropic, homoge- 
neous, Kasner cosmology H, . These initial stud- 
ies provide confidence in the technique, which is able to 
recover the homogeneity and anisotropy of the analytic 
solution, and also demonstrate both stability and second- 
order convergence to the continuum solution. 



IV. TEST-BED APPLICATIONS 

In the development of any numerical technique one 
considers a variety of test problems which provide in- 
sight into the strengths and weaknesses of the approach. 
These should be of increasing complexity, and preferably, 
cover a wide range of physical scenarios which model, in 
a simplified manner, the actual physical problems which 
motivate the code development. In this way one devel- 
ops confidence in the solutions obtained when the code 
is applied to real physical problems. Such a programme 
has been under way in the case of Regge calculus since its 
inception. Successful early applications included highly 



symmetric test problems with only a few degrees of free- 
dom, together with the first numerical construction of 
axisymmetric binary black hole initial data. The review 
by Williams and Tuckey contains a full discussion of these 
early studies j5j . 

More recently Regge calculus has been applied to in- 
creasingly complex problems of both physical and compu- 
tational import. It is these studies on which we concen- 
trate in this review. Successful applications to the Brill 
wave ^3 and black hole plus Brill wave [Tsl IT6| ini- 
tial data sets have validated computational lattice grav- 
ity, whilst also providing independent confirmation of the 
physical predictions of previous finite-difference studies. 
In this section we discuss these and other test-bed appli- 
cations of the lattice approach, allowing us to investigate 
its accuracy and flexibility. 

A. Spherical symmetry 

The Schwarzschild solution is a classic test-bed for nu- 
merical relativity. The combination of non-trivial topol- 
ogy with an event horizon covering the central singularity 
has proved a challenging problem even for modern codes. 
In this section we briefly describe an initial evolution of 
the static Schwarzschild solution using lattice gravity. 

The imposition of symmetry conditions in a lattice sim- 
ulation is challenging. While providing enormous flexi- 
bility to model complex topologies, a simplicial lattice 
(built, in three dimensions, from tetrahedra) requires 
careful construction if it is to respect the spherical sym- 
metry of the single black hole spacetime. The approach 
usually taken involves the construction of a full four- 
dimensional geometry constructed from simplices, after 
which the lattice is collapsed along the symmetry axes to 
obtain spherical symmetry in the appropriate limit. This 
approach is described elsewhere in axisymmetry '16] , and 
will also be mentioned briefly below. 




'(a) (b) 

FIG. 2: Evolution of the Schwarzschild solution in isotropic 
coordinates using a lattice of 200 vertices, (a) The radial edge 
length I (circles; every fourth point shown) plotted against 
the exponential radius coordinate rj a,t t = lOM. The exact 
solution is also shown (solid line), (b) The time evolution of 
the fractional rate of change of the radial edge lengths. After 
an initial period of relaxation, the evolution settles down to 
a static solution. 

In figure |21 we display a sample evolution of 
Schwarzschild initial data. The initial data is expressed 
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in isotropic coordinates, and the evolution is performed 
using the (fixed) analytic lapse and an area locking con- 
dition {dtgee — 0). The lattice implementation of area 
locking demands that the edges which locally span the 
0-axis remain constant, thereby fixing the area of the 
two-sphere at that vertex. After an initial relaxation 
phase, during which the maximum fractional changes in 
the edges (metric) are less than 10~^, the solution settles 
into an entirely static configuration. This application of 
lattice gravity demonstrates that the approach can yield 
accurate, stable evolutions. 

B. Axisymmetric initial data 

1. Brill waves 

The axisymmetric vacuum initial value problem first 
posed by Brill contains gravitational radiation in an 
otherwise flat initial three-geometry. Conformal decom- 
position using a flat, simply-connected base metric leads 
only to trivial solutions; Brill introduced a metric of the 
form 

dl'^^tl;^{e^'^{dp^ + dz^)+p^d(p^}, (7) 

where the arbitrary function q{p, z) can be considered 
the distribution of gravitational wave amplitude, and is 
subject to certain boundary conditions to ensure that the 
mass is asymptotically well define With this choice 
of background metric, the Haniiltonian constraint takes 
the form 

which is solved for ■0(p, z) once g(p, z) is given. 

The construction of lattice-based Brill wave initial data 
was originally considered by Dubai ^3 using a lattice 
built from prisms, and later by Gentle [T3. [la | using a 
tetrahedral lattice. The approach, modeled on the anal- 
ysis of Brill, uses a conformal decomposition technique 
and solves for the single conformal factor per vertex. The 
lattice is constructed to mirror the cylindrical polar co- 
ordinate system in which the continuum Brill wave met- 
ric is written. The conformal decomposition is obtained 
by integrating spacelike geodesies between vertices of the 
lattice, and assigning these lengths directly to the "base" 
lattice edges. 

The solutions obtained for the tetrahedral three- 
geometry were in excellent agreement with previous 
finite-difference studies Q| i thus confirming both the 
applicability of simplicial gravity, and the earlier numer- 
ical results. The prism-based calculations were consid- 
erably less accurate than the corresponding tetrahedral 
solution. Apparent horizons and ADM masses were also 
calculated for the tetrahedral initial data, again finding 
agreement with previous studies. The critical wave am- 
plitude at which an apparent horizon first forms on the 



initial surface was found to be in complete agreement 
with earlier numerical values (141 . 

2. Distorted black holes 

In this section we consider another axisymmetric, non- 
rotating configuration which contains a moment of time 
symmetry: the "distorted black hole" spacetime first con- 
sidered by Bernstein These solutions represent an 
initial slice containing a black hole together with Brill 
waves perturbations, a natural generalisation of the pre- 
ceding section, where Brill waves were considered on a 
flat Euclidean space. 

The distorted black hole spacetime is astrophysically 
interesting, as it captures the "ring-down" phase of the 
merger and coalescence of a binary black hole system. 
Following merger, a single distorted (non-spherical) black 
hole forms, then evolves towards the static Schwarzschild 
solution by the emission of gravitational radiation. This 
problem is also amenable to perturbation approaches ap- 
plied to the standard black hole solutions, providing a 
regime in which numerical relativity, theoretical analy- 
sis and gravitational wave observations combine to shed 
light on the foundations of general relativity. 

The Brill wave plus black hole spacetime UM is ob- 
tained by mirroring the original work of Brill |12| . A 
perturbation is introduced onto a "background metric", 
and the conformal factor is calculated from the single ini- 
tial value equation. Bernstein wrote the physical metric 
on the initial surface in the form 

where the exponential radius coordinate rj is defined by 
p = mexp{f])/2. The black hole topology is obtained by 
demanding that the two-sphere at = is a minimal 
surface (isometry surface), which connects two asymp- 
totically flat sheets. The "mass" m is that of the black 
hole alone, and corresponds to the ADM mass measured 
at spatial infinity when qij], 6) = 0. 

Distorted back hole initial data has been successfully 
constructed using Regge calculus 0. It is natural to 
represent the three-metric in spherical polar coordinates; 
the lattice used to model the distorted black hole space- 
time is also matched to the symmetry of the problem. 
Once again, the lattice three-geometry was obtained by 
locally aligning the background edges of the lattice with 
a spherical polar coordinate system. This implies, for 
example, that the z-axis is locally aligned with the az- 
imuthal 0-axis. This approach was used as a matter of 
convenience; there is no requirement that the lattice be 
constructed from an underlying coordinate system, but 
doing so simplifies both the application of boundary con- 
ditions and comparison with the corresponding solution 
of the Einstein equations. 

The simplicial distorted black hole initial data was 
found to agree well with finite-difference calculations, and 
estimates of the ADM mass were in excellent agreement 
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faj (b) 

FIG. 3: (a) A typical example of the lattice used to construct 
Misner initial data for a binary black hole system. The lattice 
is shown close to the z > throat in (p, z) coordinates, (b) 
The Regge conformal factor ip for the Misner binary black hole 
initial data set on the (p, z) plane with jj, — 2.2. The inner 
portion of the domain is shown. The solution is obtained on 
the first quadrant with 501 x 501 lattice vertices and an outer 
boundary at approximately — p'^ + z^ = W(f . The inset 
shows the convergence of the simplicial solution towards the 
continuum; the line of best fit has gradient 1.8. 

with previous studies. Convergence estimates indicate 
that the lattice solution converges as the second power 
of the typical lattice discretisation scale towards the true 
solution of the Einstein equations p^. 



3. Binary black holes 

Misner 17] obtained an axisymnietric solution to the 
initial value equations of general relativity suitable for 
the study of the head-on collision of equal mass, non- 
rotating black holes. Due to the complexity of the prob- 
lem there is no known analytic solution for the time de- 
velopment of this initial data. 

The maximal extension of the standard Schwarzschild 
black hole solution involves two asymptotically flat 
sheets, joined through the black hole throat. The space- 
time has an isometry through the throat. The binary 
black hole solution obtained by Misner is the nat- 
ural generalisation of this solution; it consists of two 
asymptotically flat sheets joined to one another by two 
throats. The solution represents two black holes with 
equal mass and zero angular momentum at a moment 
of time-symmetry {Kab = 0). The holes are located on 
the z-axis at 2; = ± coth /i and their throats have radius 
a — 1/ sinh/i, which also determines the ADM mass of 
each individual hole if measured in isolation. The three- 
metric takes the conformally-flat form 

df = il;^ {dp^ + dz^ + p^d^)"^) (10) 

where the conformal factor is given by 

^ ^ ( ^ 1 \ . ^ 

V' = i+V^-r^ — 7 ; + . (11) 

t{ ^1°!^ M \ + 4 ^p2 + ,2 J 



and z± = z±Zn with z„ = coth(n/u). This is a solution of 
the vacuum Hamiltonian constraint at a moment of time 
symmetry, R — 0. The single free parameter p, controls 
the separation, mass and radius of the black holes, with 
larger p, corresponding to greater separation and smaller 
bare mass [l7l| . 

In this section we present lattice solutions which cor- 
respond to the Misner initial data. Unlike the axisym- 
metric initial data presented in the previous sections, the 
existence of an analytic solution allows us to directly eval- 
uate the accuracy and convergence of the lattice solution. 
The lattice is built using Cadez coordinates which 
are spherical about both black hole throats and reduce 
to standard polar coordinates far from the origin. This 
structure necessarily introduces a coordinate singularity; 
in Cadez coordinates the singularity occurs at the origin. 

The singular nature of the "physical" coordinate sys- 
tem hindered finite-difference studies of the Misner space- 
time |l8j — we argue that the lattice approach overcomes 
this problem. While we are free to use the Cadez coor- 
dinates to construct the initial lattice ("location of ver- 
tices"; see figure ISJa)), the resulting lattice consists en- 
tirely of scalar edge lengths. Vertex positions are as- 
signed (for example) using a regularly spaced grid in 
the Cadez coordinate system, and then integration of 
geodesies is used to assign the background edge lengths 
of the lattice given the flat background metric. From 
this point on the simplicial simulation consists entirely 
of evaluating scalar functions of the edge lengths. There 
is no singularity in the lattice construction. 

A typical Regge solution for the conformal factor 
is shown in figure |3Jb), with p = 2.2. This parameter 
is used to specify the centre (2 = ±cosh/i) and radius 
(a — csch p) of the black holes. This choice guarantees 
that the initial data under consideration is of the Misner 
type, although it is clearly possible to investigate other 
configurations once the code has been benchmarked on 
the analytic solution. 

The lattice solution shown in figure Ol^b) is found to be 
in good agreement with the analytic solution, and shows 
almost second order convergence towards the analytic so- 
lution. The inset to figure Elb) shows that the conver- 
gence rate is approximately 1.8. We find that as the 
number of vertices N along each axis is increased, lattice 
edges far from the holes scale as 1/iV. Near the saddle 
point, however, the edges scale as 1/\/N. It is precisely 
this behaviour which leads to the (slightly) anomalous 
convergence rate; the slowly converging terms near the 
origin dominate. This "problem" , inherited from the sin- 
gularity in the Cadez coordinates, can be easily overcome 
by further refinement of the lattice near the saddle point. 

Despite these minor issues, the Regge approach recov- 
ers a second-order accurate approximation to the solution 
of Einstein's equations for the Misner binary black hole 
initial data. Following further refinement of the lattice 
near the saddle point, the method should yield precisely 
second order convergence. The Cadez coordinates are 
used only as a convenience when constructing the initial 
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FIG. 4: Brill wave evolution, (a) Time evolution oi K — K\, 
evaluated at the origin. We use the fact that K is equivalent 
to the fractional rate of decrease of the volume per vertex. 
Evolution is shown for a wave of amplitude a = 0.01 on a 
100 X 100 grid, timestep At — Ap/4, and with the outer 
boundary at p = 2; = 10. (b) Mean fractional difference 
between magnitude of the first "bump" in (a) and the cor- 
responding "continuum value", plotted as a function of the 
grid resolution N x N. The continuum value is estimated us- 
ing Richardson extrapolation of the two most accurate values. 
We observe close to second order convergence; the gradient of 
the hue of best fit is approximately 1.9. 

lattice. If and when such an approach proves problem- 
atic, the lattice can, and should, be further subdivided 
to obtain a consistent "triangulation" . 

C. Brill wave evolutions 

In this section we describe an application of lattice 
gravity to the evolution of the (low amplitude) Brill wave 
initial data discussed above. Although all calculations 
are performed using lattice gravity, in this section we de- 
scribe the techniques and gauge choices in the standard 
language of numerical relativity. Although the Regge 
framework is complete and self-contained, we find it use- 
ful to highlight connections with the continuum and more 
"standard" approaches to numerical relativity. 

The spatial metric is written in cylindrical polar coor- 
dinates, and assumed to take the form 

ds^ = V'^ { (V + dz^) + p'^dcj)'^ } (12) 

throughout the evolution, where ^ = ^{p,z,t) and -0 — 
tp{p,z,t). The lattice is adapted to these coordinates, as 
described elsewhere This form of the three-metric 

7ab was also chosen by Garfinkle and Duncan and 
implicitly involves the imposition of the two conditions 

7pp = Izz, and = 0, (13) 

in addition to the initial assumption of axisymmetry. 
These two conditions are effectively a choice of gauge, 
and through the Einstein equations, determine the shift 
vector /3* [13. Whilst the same holds true in the lattice 
approach, we follow a different strategy in this initial 
evolution of gravitational radiation on a lattice. 

The conditions (|13|l are applied directly to the lat- 
tice edges corresponding to the spatial components of 



the three-metric. A zero shift condition, /3* = 0, is also 
implemented on the lattice to ensure that the timelike 
worldlines generated by the vertices are orthogonal to the 
lower hypersurface. This is an over-specification of the 
available degrees of freedom, but we argue that for weak 
Brill waves (a << 1) the spacetime is a perturbation of 
flat space, and thus all of the above conditions can be 
enforced to within a reasonable accuracy. The choice of 
zero shift, together with the extra gauge choice H13|l . will 
in general violate the Einstein equations and their simpli- 
cial counterparts. For small perturbations away from flat 
space we expect that these errors are controllable. The 
generic evolution of medium to large amplitude waves will 
require a consistent specification of the available gauge 
freedoms. The choices outlined above, however, are suffi- 
cient to perform the initial low-amplitude test evolutions 
presented in this paper. 



The remaining gauge freedom is determined by setting 
a = 1 (unit lapse). This condition is not ideal, as it is 
known that even weak waves will impart non-zero veloc- 
ities to the vertices of the lattice (or grid). Once the 
wave has propagated away from the centre of the grid we 
expect to recover a non-trivial, non-static coordinatisa- 
tion of flat space. For strong amplitude waves we will 
require more advanced slicing conditions, such as max- 
imal (Tri^T = 0) or harmonic slicing. Geodesic slicing 
is sufficient, however, for this initial low-amplitude test 
problem. 



Boundary conditions are imposed on the lattice to en- 
force the reflection symmetry about the p — Q axis, to- 
gether with reflection symmetry about the z = axis to 
simplify the computations. These conditions are imple- 
mented using ghost cells centred about the axis. This 
approach has been used successfully in previous studies 
[l9j , and is found to work well. The outer boundary con- 
ditions arc implemented using a radiative "Sommerfeld" 
condition, implemented in differential form on and ^ 
along the outer boundary • 

Results are shown in figures 01 and |S1 and demonstrate 
that Regge calculus can successfully and accurately prop- 
agate a gravitational wave through the lattice. An initial 
convergence estimate is performed by examining the con- 
vergence of the magnitude of the first "bump" in figure 
0Ja). This feature of the solution is found to converge at 
close to second order, demonstrating a degree of consis- 
tency in the solutions obtained. Figure [3 displays the de- 
viation of the conformal factor [if)'^ — 1) at various times; 
the propagation of the wave can be observed, with the 
inner region returning to a roughly constant value ( "flat- 
ness" ) once the wave has escaped towards the boundary. 
Reflections from the outer boundary currently prevent 
long term evolutions. 
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FIG. 5: The deviation of the conformal factor, tp^ — 1, plotted 
as a function of p along the line z = for various times {(a) 
t = to 0.5, (h) t = 1 to 2, fcj t = 3 to 4. fdj t = 5 to 
6) for a Brill wave with amplitude a = 0.01. We use the 
"Eppley-type" Brill function q{p,z) |20|. and show solutions 
for a lattice with 100 x 100 vertices. Note that the x-axis 
scale changes in parts (c) and (d), and the outer boundary is 
placed at p = 2 = 10. 



V. THE FUTURE OF LATTICE APPROACHES 
TO NUMERICAL RELATIVITY 

We have reviewed several recent applications of Regge 
calculus, and demonstrated that it provides a unique, 
and thus far successful, alternative to the more standard 



techniques employed in numerical relativity. 

In particular, we have presented new solutions of lat- 
tice gravity corresponding to initial data for the head-on 
collision of black holes, together with the first successful 
time evolution of gravitational radiation on a lattice. De- 
spite these and the other studies described above, signif- 
icant questions remain regarding the fundamental struc- 
ture of Regge calculus. 

Much work remains to be done to address the dual 
and diffeomorphic structures of the lattice, the relation 
of this and other lattice approaches to standard finite 
element and finite volume discretisations of differential 
equations, and the inclusion of matter. On the numerical 
front, further applications in three-plus-one dimensions 
are vital, requiring the development of a modern parallel 
code. 

Despite these issues, every indication is that the lat- 
tice approached to gravity developed by T. Regge will 
continue to provide an interesting and complementary 
approach to numerical relativity. 
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